First-order chiral transition in the compact lattice theory of graphene 
and the case for improved actions 
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A comparison of the compact and non-compact lattice versions of the low-energy theory of 
graphene is presented. The compact theory is found to exhibit a chiral phase transition which 
appears to be of first order, at a critical coupling of fi c = 0.42 ± 0.01. We confirm that the non- 
compact theory exhibits a second-order transition at /3 C = 0.072 ±0.003, and determine the effects of 
UV-divergent tadpole contributions in both cases. Upon tadpole improvement of the non-compact 
theory we find fjj 1 = 0.163 ± 0.002, which strengthens the case for a semimetal- insulator transition 
in graphene at strong Coulomb coupling. Finally, we highlight the need for systematic studies using 
improved lattice actions. 

PACS numbers: 73.63.Bd, 71.30.+h, 05.10.Ln 



Graphene, a sheet of sp 2 -bonded carbon, has become 
an attractive candidate for nanoscale electronics due 
to its many remarkable and desirable properties [l|, Q. 
These include high carrier mobility at room temper- 
ature, great mechanical and tensile strength, in addi- 
tion to chemical stability and impermeability. While the 
hexagonal lattice symmetry suggests that graphene is a 
semimetallic material exhibiting massless Dirac quasipar- 
ticles, the possibility that graphene might become an ex- 
citonic insulator due to strong Coulomb interactions has 
recently been revived 

In a series of papers within the Lattice Monte 
Carlo (LMC) framework @, the semimetal-insulator 
transition, which manifests itself as the spontaneous 
breaking of a U (4) chiral symmetry, was found to happen 
at a critical coupling of a c ~ 1.1, which is intermediate 
between that of graphene on a Si0 2 substrate (a ~ 0.8) 
and suspended graphene (a ~ 2.1). LMC studies of 
a closely related Thirring-like model, including the de- 
termination of the renormalized Fermi velocity, were re- 
ported in Ref . Q . These results belong to a larger class of 
LMC studies, such as those of Quantum Electrodynam- 
ics in 2 + 1 (QED 3 ) and 3 + 1 (QED 4 ) H dimensions, 
and of four-fcrmion theories such as the Thirring [9j and 
Gross-Neveu [l(| models. 

Spontaneous chiral symmetry breaking in LMC is typi- 
cally studied using staggered fermions 11], as chiral sym- 
metry is then partially preserved at finite lattice spacing 
a. In practical LMC simulations the continuum limit 
is recovered in the vicinity of second-order phase tran- 
sitions, where the relevant correlation lengths £ diverge. 
In the case of graphene, such a continuum description is 
attainable on the critical line of vanishing bare mass and 
(inverse) coupling /3 > /3 C ~ 0.072, with a c = 1/(4tt/3 c ). 
In the strong coupling phase /3 < /3 C the Coulomb in- 
teraction induces the formation of particle-hole pairs, 
with a binding energy that yields an intrinsic cutoff scale 



rather, it defines an unstable theory). Keeping a finite, 
however, one may approach the critical line and recover a 
well-defined theory, with bound states at strong coupling 
and massless fermions at weak coupling. 

The requirement of exact gauge invariance on the lat- 
tice necessitates the use of "gauge links" U = exp(i6), 
with 9 the lattice gauge field. Gauge links have the side 
effect of introducing vertices of higher order in a, such 
as photon-photon interactions, which are absent in the 
continuum theory. Such vertices also yield potentially 
large "tadpole" contributions, where the naive power- 
counting in a is cancelled by UV divergences. Conse- 
quently, Lattice QCD simulations employ various levels 
of "improvement" 12J, in order to minimize the impact 
of such discretization artifacts. We report here a sys- 
tematic study of the tadpole effects in the graphene the- 
ory, which were first pointed out in Ref. 13], and their 
impact on the determination of (3 C . We also present a 
comparison of the non-compact and compact lattice the- 
ories of graphene, where the latter case involves photon 
self-interactions which are absent in the former. 

The Euclidean action of the lattice theory of graphene 
is conventionally split into gauge and fermion compo- 
nents, such that S E = S 9 E + S E . In the compact formu- 
lation, the gauge field 9 enters into S E in terms of link 
variables, giving 
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such that the limit a — ¥ is not well defined (or 



where K(x) denotes the real part of x, (3 = l/g 2 is the 
inverse coupling, n a site on the (3+l)-dimensional space- 
time lattice, and a unit vector in the direction /i. 

Early on, the analogous compact formulation of QED 4 
was found to exhibit a first-order transition [l4| where 
the breaking of chiral symmetry is coincident with the 
condensation of magnetic monopoles. While the addition 
of a four-Fermi interaction to QED 4 makes it possible 
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to isolate the chiral symmetry breaking transition, the 
resulting theory may belon g, to a different universality 
class than continuum QED [l5|. In order to avoid this 
situation for the case of graphene, an attractive option is 
the non-compact formulation 
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which is free of photon self-interactions, and is known to 
support a second-order transition which may be identi- 
fied with the spontaneous breaking of chiral symmetry 
in the continuum theory. Thus, Eq. ^ has become the 
standard choice for LMC studies of Abelian gauge the- 
ories, such as QED 3 Q, QED 4 Q and we have used it 
for graphene. The first objective of this study is to char- 
acterize the compact version of the low-energy theory of 
graphene, which differs significantly from QED 4 as the 
spatial gauge links are constant and the fermions propa- 
gate in two spatial dimensions only. 

The fermion action is identical in the compact and 
non-compact theories, and in the staggered fermion for- 
mulation it is given by 



(3) 



where (n, n') denote the sites of a (2+l)-dimensional 
space-time sublattice, and gauge invariance is retained 
by coupling the staggered spinors x n to the gauge field 
via link variables in the time direction. The staggered 
Dirac operator is given by (see e.g. Ref. flij ]) 
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where rj 1 n = (— l) n ° and rj 2 n = (— l)"°+ n i. The Fermi 
velocity v can be absorbed into the remaining param- 
eters, giving /3 = v/g 2 and m = m /v. Simulations 
are thus conventionally performed in terms of (/3, m) and 
v = 1. The mass term acts as a symmetry breaking pa- 
rameter, without which the chiral condensate a would 
be zero at finite volume. The limit m — > is reached 
by extrapolation, using an equation of state (EOS) of 
the form m = f(a,0), which describes a second-order 
transition with critical exponents 5 and (3. A detailed 
description of the EOS can be found in Ref. [5[ . 

"Tadpole improvement" (TI) is a non-perturbative 
method due to Lepage and Mackenzie [17| that accounts 
for the UV divergent tadpole contributions by renormal- 
izing the link field U. The effect of the UV modes is 
encoded in the function u , which depends on the input 
parameters of the simulation and is to be determined a 
posteriori. A conventional definition of u is 

U Q = {P) X '\ P=^J2 U n U t+^ ( 5 ) 
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FIG. 1: (Color online) Plaquettes (P) as a function of j3 
for the non-compact (left) and compact (right) theories, for 
Eq. ([5J averaged over the fermionic sublattice. Our notation 
is N z x N% , where N x is the extent of the fermion sublattice 
and N z that of the bulk dimension. Black datapoints are for 
a 32 4 lattice, while red (gray) points are for 12x28 3 . To illus- 
trate the depencence on m, the datapoints have been shifted 
horizontally relative to m = 0.010, and for the compact case 
the 24 4 data have also been shifted vertically relative to the 
32 4 data. 



in terms of the plaquette P. It should be noted that the 
power 1/2 (instead of 1/4 as in Lattice QCD) is due to 
the smaller number of fluctuating gauge links. 

Recently, Ref. [l3| has reported that it (/3, m) deviates 
significantly from unity, and thus the effects of TI are 
likely to be significant. Therefore, the second objective 
of this study is to determine how TI affects the deter- 
mination of (5 C and the critical exponents. In the non- 
compact theory, modifications due to TI are restricted 
to the fermion action, whereas in the compact case the 
gauge action is affected as well. 

As in Ref. (l3j . we consider the plaquettes in the 
(x, t) and (y, t) planes, however unlike Ref. [l3[ we com- 
pute (P) by summing over the plaquettes in the (2+1)- 
dimensional fermionic sublattice only. This choice is ap- 
propriate for the non-compact theory, and the resulting 
numerical differences arc insignificant for the subsequent 
analysis. Our results for P are shown in Fig. [TJ We find 
that P is independent of m up to statistical fluctuations, 
and thus we define u (/3) as the average of P over m. 

In practice, TI amounts to replacing the link field U 
in Eqs. JT]) and ([3]) according to U — > U/u , where u 
is to be determined a posteriori using Eq. ([5]). If the 
staggered spinors are rescaled as x E 
define an "improved" Dirac operator 
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FIG. 2: (Color online) Chiral condensate a (left) and sus- 
ceptibility Xi = da /dm (right) in the non-compact theory, 
together with an EOS fit to the unimproved lattice data. The 
fitted parameters are /3 C = 0.072 ± 0.003 and 8 = 2.3 ± 0.3, 
with f3 ~ 1. The data are obtained for lattices of size 12 x 32 3 
and 32 4 (circles) and 12 x 28 4 (squares). Errors for the in- 
dividual datapoints were obtained with the standard block- 
jackknife method. Note the deviations from the EOS (scaling 
violations) at small ft and large m. 



FIG. 3: (Color online) Chiral condensate a (left) and suscep- 
tiblity xi (right) in the non-compact theory after TI, together 
with an EOS fit. Note the significantly reduced scaling vio- 
lations compared to the unimproved data. The optimal fit 
is compatible with = 1, giving = 0.163 ± 0.002 and 
S = 2.2 ± 0.1. The errors are mostly systematical, due to 
finite- volume effects and residual scaling violations. The black 
symbols denote data on 12 x 32 3 and 32 4 lattices, while red 
(gray) symbols denote data for 12 x 28 3 . 



such that the tadpole-improved observables may be com- 
puted using the original unimproved gauge configura- 
tions, provided that the input parameters of the sim- 
ulation and the spinors are reinterpreted according to 
Eq. J7]). As m < 1, the net effect of TI (apart from 
possible shifts of the critical coupling f3 c ) is to make the 
condensate smaller in the spontaneously broken phase, an 
effect which increases with decreasing f3. For the compact 
gauge action in Eq. ([1} , TI leads to a similar prescription 
g' = u <7, while the non-compact action in Eq. ^ is not 
directly modified as it involves no gauge links. 

The primed parameters in Eq. ([7]) represent the input 
for the LMC calculation. Once u (f3) has been mapped 
out, the tadpole-improved (unprimed) quantities can be 
determined. Apart from the rescaling of the chiral con- 
densate according to a' u a, we find for the compact 
case 



v'/u 



g 2 g' 2 /ul 
while for the non-compact case we have 
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and we note that m = m /v remains unmodified. Both of 
these results differ from the prescription ft = u 2 /3' which 
was applied in Ref. [13]. However, such a choice would be 
valid for the compact theory in the absence of TI for the 
fermion action. It is noteworthy that the non-compact 
theory would have j3' = (3 if the tadpoles of fermionic 



origin were ignored. Our results for the non-compact 
theory are shown in Figs. [2] and [3] for the unimproved 
and TI cases, respectively. Similarly, the results for the 
compact theory are given in Figs. 0] and [3] Most of our 
data for the non-compact case are taken from Ref. [j|, 
except for the data on 32 4 lattices close to the critical 
point. 
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FIG. 4: (Color online) Chiral condensate a (left) and suscep- 
tibility Xi (right) for the compact, unimproved theory as a 
function of /3. A very sharp, possibly first-order transition is 
observed at /3 C = 0.42 ±0.01. Data on a 24 4 lattice are shown 
in red (gray) and are for clarity slightly offset from the (black) 
datapoints on a 32 4 lattice. The volume dependence of Xi is 
greatly enhanced in the spontaneously broken phase. 
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FIG. 5: (Color online) Chiral condensate a and susceptibility 
Xi in the compact theory of graphene after TI. The presenta- 
tion is similar to Fig.0] The critical coupling ^J 1 = 0.27±0.01 
has been shifted closer to that of the non-compact theory. At 
P c , the data indicate a discontinuous jump of A CT ~ 0.04 in 
the condensate. 



As the tadpole correction u (f3) ~ 0.5 close to the chi- 
ral phase transition in the non-compact theory, the effects 
of TI can potentially be dramatic. We find that the agree- 
ment with the EOS improves and the effects of outlying 
data points are lessened, which is due to decreased scaling 
violations at small (3. While the net effect on the critical 
exponents is slight (these remain compatible with the val- 
ues of Ref. |5j), the change in /3 C is larger as a consequence 
of Eq. ©. After TI, we find /3J 1 ~ 0.16, which is a factor 
of ~ 2 larger than the unimproved value. On the other 
hand, the compact theory exhibits an abrupt transition 
close to /3 ~ 0.45, and is furthermore not compatible 
with the EOS description of Ref. @, although TI does 
bring j3 c closer to the non-compact result, as evidenced 
by Fig. [SJ As a first estimate, we find A CT ~ 0.04 ± 0.01 
for the jump in the condensate. 

It is intriguing that /3J 1 > f3 c in the non-compact the- 
ory, as it suggests that graphene samples on a Si0 2 sub- 
strate may become insulating as well. However, LMC 
studies of the Fermi velocity in graphene @ indicate a 
significant downward renormalization of v due to strong 
Coulomb interactions, thereby again decreasing the phys- 
ical value of j3 c . Further LMC studies of the electrical 
conduction properties of graphene are in progress (lij . 

In conclusion, we find that the effect of TI is not as 



drastic as reported in Ref. 13], particulary in the non- 
compact theory, where f3 c is shifted to a larger (instead 
of a smaller) value. The predicted semimetal-insulator 
transition in suspended graphene is therefore unlikely to 
be an artifact of tadpole effects. While the full chiral 
symmetry is not realized in LMC simulations with stag- 
gered fermions at small /? (see Refs. [1, EH), it is known 
to eventually be restored in the vicinity of the tricritical 



ness may be inconsequential, as the accuracy of the ex- 
trapolation can be systematically improved by obtaining 
additional data closer to the tricritical point. Neverthe- 
less, simulations with overlap fermions [19j may serve to 
further clarify this point, and to extend the scope of the 
present work. 
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